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Abstract 

In this paper we construct spherical thin-shell wormholes supported by a Chaplygin gas. 
For a rather general class of geometries we introduce a new approach for the stability analy- 
sis of static solutions under perturbations preserving the symmetry. We apply this to worm- 
holes constructed from Schwarzschild, Schwarzschild-de Sitter, Schwarzschild-anti de Sitter and 
Reissner-Nordstrom metrics. In the last two cases, we find that there are values of the param- 
eters for which stable static solutions exist. 
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^ ■ 1 Introduction 

m 

Traversable Lorentzian wormholes [1] are solutions of the equations of gravitation associated to 
a nontrivial topology of the spacetime: their basic feature is that they connect two regions (of 
the same universe or two separate universes [1,2]) by a throat. For the case of static wormholes 
the throat is defined as a minimal area surface satisfying a flare-out condition [3]. To fulfill this, 
wormholes must be threaded by exotic matter that violates the null energy condition [1-4]; it was 
shown by Visser et al [5], however, that the amount of exotic matter needed around the throat can 
be made as small as desired by means of an appropriate choice of the geometry of the wormhole. 



A well studied class of wormholes is that of thin-shell ones, which are constructed by cutting 
and pasting two manifolds [2,6] to form a geodesically complete new one with a shell placed in the 
joining surface. This makes such wormholes of particular interest because the exotic matter needed 
for the existence of the configuration is located only at the shell. Stability analysis of thin-shell 
wormholes under perturbations preserving the original symmetries has been widely developed. A 
linearized analysis of a thin-shell wormhole made by joining two Schwarzschild geometries was 
performed by Poisson and Visser in Ref. [7]. Later, the same method was applied to wormholes 
constructed using branes with negative tensions in Ref. [8] , and the case of transparent spherically 
symmetric thin-shells and wormholes was studied in Ref. [9]. The linearized stability analysis 
was extended to Reissner-Nordstrom thin-shell geometries in Ref. [10], and to wormholes with a 
cosmological constant in Ref. [11]. The case of dynamical thin-shell wormholes was considered in 
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Ref. [12]. The stability and energy conditions for five dimensional thin-shell wormholes in Einstein- 
Maxwell theory with a Gauss-Bonnet term were studied in Ref. [13], while thin-shell wormholes 
associated with cosmic strings have been treated in Refs. [14]. Other related works can be found 
in Refs. [15]. 

The requirement of matter violating energy conditions relates the study of wormholes to mod- 
ern cosmology: Current day observational data seem to point towards an accelerated expansion of 
the universe [16]. If General Relativity is assumed as the right gravity theory describing the large 
scale behavior of the universe, this implies that its energy density p and pressure p should violate 
the strong energy condition. Several models for the matter leading to such situation have been 
proposed [17]. One of them is the Chaplygin gas [18], a perfect fluid fulfilling the equation of state 
pp = —A, where A is a positive constant. A remarkable property of the Chaplygin gas is that the 
squared sound velocity v 2 = A/p 2 is always positive even in the case of exotic matter. Though 
introduced for purely phenomenological reasons (in fact, not related with cosmology [19]), such an 
equation of state has the interesting feature of being derivable from string theory; more precisely, it 
can be obtained from the Nambu-Goto action for d— branes moving in {d + 2)— dimensional space- 
time if one works in the light-cone parametrization [20]. Besides, an analogous equation of state, 
but with A a negative constant, was introduced for describing cosmic strings with small structure 
("wiggly" strings) [21]. 

Models of exotic matter of interest in cosmology have already been considered in wormhole con- 
struction. Wormholes supported by "phantom energy" (with equation of state p = up, u < — 1) 
have been studied in detail [22]. A generalized Chaplygin gas, with equation of state pp a = —A 
(0 < a < 1), has been proposed by Lobo in Ref. [23] as the exotic matter supporting a wormhole 
of the Morris-Thorne type [1]; there, as a possible way to keep the exotic matter within a finite 
region of space, matching the wormhole metric to an exterior vacuum metric was proposed. If, 
instead, a thin-shell wormhole is constructed, exotic matter can be restricted from the beginning 
to the shell located at the joining surface. In the present paper we study spherically symmetric 
thin-shell wormholes with matter in the form of the Chaplygin gas (the generalized Claplygin gas 
introduces a new constant and non trivial complications in the equations, which possibly cannot 
be solved in an analytical way). We introduce a new approach for the study of the stability under 
radial perturbations. The more complex case of the stability analysis of thin-shell wormholes under 
perturbations that do not preserve the symmetry has not been addressed in previous works, even 
for simpler metrics and equations of state, so we consider it beyond the scope of this article. In 
Section [2] we apply the Darmois-Israel formalism to the cut and paste construction of a generic 
wormhole with the Chaplygin equation of state imposed on the matter of the shell. In Section Owe 
perform a detailed analysis of the stability under spherically symmetric perturbations. In Section 
d]we analyze the specific cases of the Schwarzschild and the Reissner-Nordstrom geometries, and 
we also consider the inclusion of a cosmological constant of arbitrary sign. In Section [5] the results 
are discussed. We adopt units such that c = G = 1. 



2 Wormhole construction 

Let us consider a spherically symmetric metric of the form 

ds 2 = -f{r)dt 2 + /{r^dr 2 + r 2 {d9 2 + sin 2 Odyf 2 ), (1) 
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where r > is the radial coordinate, < 6 < tt and < <p < 2ir are the angular coordinates, and 
f(r) is a positive function from a given radius. For the construction of the thin-shell wormholes, 
we choose a radius a, take two identical copies of the region with r > a: 

M ± = {X a = (t,r,0^)/r>a}, (2) 

and paste them at the hypersurface 

S = S ± = {X/F(r) = r - a = 0}, (3) 

to create a new manifold M = M + U M~ . If the metric ([T]) has an event horizon with radius 
r/,, the value of a should be greater than r^, to avoid the presence of horizons and singularities. 
This construction produces a geodesically complete manifold, which has two regions connected by 
a throat with radius a, where the surface of minimal area is located and the condition of flare- 
out is satisfied. On this manifold we can define a new radial coordinate I = db J \J\j f(r)dr 
representing the proper radial distance to the throat, which is situated at / = 0; the plus and 
minus signs correspond, respectively, to Ai + and Mr . We follow the standard Darmois-Israel 
formalism [24, 25] for its study and we let the throat radius a be a function of time. The wormhole 
throat X is a synchronous timelike hypersurface, where we define coordinates ^* = (r, 6, ip), with r 
the proper time on the shell. The second fundamental forms (extrinsic curvature) associated with 
the two sides of the shell are: 



(4) 



where are the unit normals (n 7 re 7 = 1) to S in A4: 



3 a dF OF 



-1/2 Qp 

Wr- (5) 



Working in the orthonormal basis {ef,eg,e^} (ef = e T , eg = a 1 eg, = (asin^) 1 e ¥ ,), for the 
metric (pQ) we have that 

K% = K% = ±iy/f(a) + a*, (6) 



1 

a 

and 



± _ f(a) + 2a 

2 v/( a ) + a 

where a prime and the dot stand for the derivatives with respect to r and r, respectively. Defining 
[K,.] = Kt — Kz, K = tr[K$j\ = [K l ~] and introducing the surface stress-energy tensor S r = 
diag(a, Pg,p0) we obtain the Einstein equations on the shell (the Lanczos equations): 

- [K i} ] + K m = 8vr%, (8) 

which in our case correspond to a shell of radius a with energy density a and transverse pressure 
P = Pg=P<P given by 

V/(a)+« 2 , (9) 



a 



2ira 



1 2aa + 2d 2 + 2/(q) + af'(a) _ 

p= = . 10 

8vra V/(a) + a 2 
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The equation of state for a Chaplygin gas has the form 

-A , , 

P= — , (11) 
a 

where A is a positive constant. Replacing Eqs. Q and (Tl"0j) in Eq. (fTTI) . we obtain 

2ad + 2d 2 - 16vr 2 Aa 2 + 2/(a) + a/'(a) = 0. (12) 

This is the differential equation that should be satisfied by the throat radius of thin-shell wormholes 
threaded by exotic matter with the equation of state of a Chaplygin gas. 

3 Stability of static solutions 

From Eq. (|12p . the static solutions, if they exist, have a throat radius ao that should fulfill the 
equation 

- 16vr 2 Aa 2 Q + 2/(a ) + a /'(a ) = 0, (13) 

with the condition ao > r h if the original metric has an event horizon. The surface energy density 
and pressure are given in the static case by 



and 

2irAa 
V/( a o) 

The existence of static solutions depends on the explicit form of the function /. To study the 
stability of the static solutions under perturbations preserving the symmetry it is convenient to 
rewrite a(r) in the form 

a(r) = ao [l + e(r)], (16) 
with e(r) <C 1 a small perturbation. Replacing Eq. (fl"6|) in Eq. (JT2"j) and using Eq. (fl3|) . we obtain 

(1 + e)e + e 2 - %k 2 A{2 + e)e + g{a , e) = 0, (17) 

where the function g is defined by 



2f(a + a e) + a (l + e)/'(a + a e) 2f(a ) + a f'(a ) 
5(ao ' £) = 2^ 2~4 • 

Defining z/(r) = e(r), Eq. ()17p can be written as a set of first order differential equations 

v 

. _ 8-k 2 A(2 + e)e - g(a ,e) v 

TT~e TT7 

Taylor expanding to first order in e and v we have 



e = v 

2 AfO I r,lr, ^\ ,,2 

(19) 



e = v 

v = Ae, (20) 
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where 

A = l^A - & l (a Q , 0) = l^A - 3/ ' (ao) + ao///(Qo) , (21) 

which, by defining 

) andM=f ° M, (22) 



^ / V A , 

can be put in the matrix form 

i = M£. (23) 

If A > the matrix M has two real eigenvalues: Ai = — vA < and A2 = vA > 0. The presence 
of an eigenvalue with positive real part makes this case unstable. As the imaginary parts of the 
eigenvalues are zero, the instability is of saddle type. When A = 0, we have Ai = A2 = 0, and to 
first order in e and v we obtain v = constant = vq and e = eo + vq (r — To), so the static solution is 
unstable. If A < there are two imaginary eigenvalues Ai = — iv/jAj and A2 = i^/|A|; in this case 
the linear system does not determine the stability and the set of nonlinear differential equations 
should be taken into account. For the analysis of the A < case, we can rewrite Eq. (|19p in polar 
coordinates (p, 7), with e = pcos7 and v = psin7, and make a first order Taylor expansion in p, 
which gives 

p = shi7cos7(l + A)p 

7 = A cos 2 7 — sin 2 7 + h(j)p, (24) 

where ^(7) is a bounded periodic function of 7. For small values of p, i.e. close to the equilibrium 
point, the time derivative of the angle 7 is negative (the leading term A cos 2 7 — sin 2 7 is negative), 
then 7 is a monotonous decreasing function of time, so the solution curves rotate clockwise around 
the equilibrium point. To see that these solution curves are closed orbits for small p, we take a time 
T\ so that (e(n), = (ei,0) with e\ > 0. As the solution curve passing through (ei,0) rotates 

clockwise around (0, 0), there will be a time T2 > t\ such that the curve will cross the e axis again, in 
the point (e(r2), v{j<if) = (£2, 0), with 62 < 0. As the Eq. (fT9l) is invariant under the transformation 
composed of a time inversion r — > —t and the inversion v — * —v, the counterclockwise curve 
beginning in (ei,0) should cross the e axis also in (e2,0). Therefore, for A < the solution curves 
of Eq. (fT9|) should be closed orbits near the equilibrium point (0,0), which is a stable center. The 
only stable static solutions with throat radius ao are then those which have A < 0, and they are 
not asymptotically stable, i.e. when perturbed the throat radius oscillates periodically around the 
equilibrium radius, without settling down again. 

4 Application to different geometries 

In this section we analyze wormholes constructed with different metrics with the form of Eq. ((T|) . 

4.1 Schwarzschild case 

For the Schwarzschild metric we have that 

r 

where M is the mass. This geometry has an event horizon situated at = 2M, so the radius of 
the wormhole throat should be taken greater than 2M. The surface energy density and pressure 
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Figure 1: Solutions of Eq. (f2"8|) as functions of AM 2 . For AM 2 < (547r 2 ) -1 there are three real 
roots, two of them positive Oq (full black line) and (full grey line), and one negative a^ 9 (dashed 
line); when AM 2 = (54-7r 2 )~ 1 the two positive roots merge into a double one; and if AM 2 > (547r 2 )~ 1 
there is only one negative real root. 



for static solutions are given by 



G = 3/2 > ( 26 ) 

Z7ra 



and 

P= 2 ' A "° 2 , (27) 

with the throat radius ao that should satisfy the cubic equation: 

8Tr 2 Aal - a + M = 0. (28) 

Using that A > and M > 0, it is not difficult to see that this equation has one negative real root 
and two non real roots if AM 2 > (547r 2 ) -1 , one negative real root and one double positive real root 
if ^4M 2 = (54-7T 2 ) -1 , and one negative real root and two positive real roots if AM 2 < (54-7r 2 )~ 1 . 
The solutions of cubic equation ((2*51) , plotted in Fig. [TJ are given b}0 



neg 



2/3 

-1 - iy/3 - (1 - iy/3) ( -3ttV6AM + iy/l - 54tt 2 AM 2 " 



4it\/6A (-3it\/6AM + iy/1 - 5An 2 AM 2 ^j 



1/3 



(29) 



1 + iy/3 - (l + iy/3) SttVMM + iy/1 - 54ir 2 AM 2 



2/3 



< = * -m — '— , (30) 

Airy/M ( -SttVMM + iy/1- Mir 2 AM 2 



1 the notation for the roots will be clear below 
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and 

/ , , \ 2/3 

1 + -SttVgAM + Wl - 54it 2 AM 2 

°0 = ^7 — \ l/3 > (::! " 



where the powers of complex numbers give the principal value. The negative root a^ 9 has no 
physical meaning, so if AM 2 > (547r 2 ) -1 there are no static solutions. When AM 2 = (54-7T 2 ) -1 the 
positive double real solution of Eq. (f28j) is clq = 3M/2 < r^, then no static solutions are present. 
For AM 2 < (547T 2 )- 1 , the positive roots of Eq. (|28|) are M < a s < 3M/2 and 3M/2 < a%; the first 
one is always smaller than r^, thus it has to be discarded, and the second one is greater than if 
AM 2 < ao, where ao ~ 1.583 x 1CP 3 < (547T 2 ) -1 (obtained numerically). To study the stability of 



Oq, we calculate A for this case: 



A = 16vr 2 ^ 5-, 32 

a* 



which, with the help of Eq. ([28jl . can be simplified to give: 

A = ^r(2a -3M). (33) 
% 

Then, using that Oq > 3M/2, it is easy to see that A is always a positive number so, following Sec. 
[3l the static solution is unstable (saddle equilibrium point). Briefly, for the Schwarzschild metric 
if AM 2 > ao no static solutions are present, and if AM 2 < uq there is only one unstable static 
solution with throat radius clq. 

4.2 Schwarzschild-de Sitter case 

For the Schwarzschild-de Sitter metric the function / has the form 

/(,-, V, (34) 

r 3 

where A > is the cosmological constant. If AM 2 > 1/9 we have that f(r) is always negative, so we 
take < AM 2 < 1/9. In this case the geometry has two horizons, the event and the cosmological 
ones, which are placed, respectively, at 



1 + *V3 - (1 + *V3) (sVAM + iy/l - 9AM 2 ) 



2/3 

rh = " i/ 3 " ' ( 35 ) 

2yfA (-3VAM + iy/l - 9AM 2 

/ ^_ , \ 2/3 

1 + -3VAM + Wl - 9AM 2 

r c = ) (36) 

\/A (sVAM + Wl - 9AM 2 

The event horizon radius is a continuous and increasing function of A, with r^(A — > + ) = 2M 
and rh(AM 2 = 1/9) = 3M, and the cosmological horizon radius r c is a continuous and decreasing 
function of A, with r c (A -> 0+) ->• +oo and r c (AM 2 = 1/9) = 3M. If < AM 2 < 1/9 the wormhole 
throat radius should be taken in the range r/j < ao < r c , and if AM 2 = 1/9 the construction of the 
wormhole is not possible, because = r c = 3M. Using Eqs. (|14p and (|15p . we obtain that the 
energy density and the pressure at the throat are given by 



V-Aog + 3a - 6M 



^=- v ; r , 3 2 — . ( 37 ) 
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Figure 2: Solutions of Eq. (|39|) as functions of AM 2 , where A = A + A/(127r 2 ). We have that A 
is always positive or it has any sign, depending on A > or A < 0, respectively. When A < 0, 
there is only one real (and positive) root Oq (full grey line); for < AM 2 < (547T 2 ) -1 there are 
three real roots, two of them positive Oq (full black line) and (full grey line), and one negative 
dp 69 (dashed line); when AM 2 = (547r 2 ) -1 the two positive roots merge into a double one; and if 
AM 2 > (547r 2 )~ 1 there is only one negative real root. 



2\/3ttA 



3/2 



p= . V u (38) 
y'-Aog + 3a - 6M 



The throat radius ao should satisfy in this case the cubic equation 

8^ 2 \ A + ^2 ) a o - °o + M = 0. (39) 

If we define A = A + A/(127r 2 ), which in this case is a positive number, it is easy to see that Eq. 
([39]) has the same form as Eq. (j28j) . so the solutions of Eq. ([39]) are given again by Eqs. (|29M3ip . 
with A replaced by A. These solutions are shown in Fig. [2] (the part of the plot with A > 0). 
With the same arguments of Sec. 14.11 if AM 2 > (547r 2 )~ 1 we have no static solutions. Using that 
2M < r/j < 3M, when AM 2 < (54-7T 2 ) -1 , aff is always smaller than rh, so it has to be discarded, 
and Oq can be greater, equal or smaller than rh, depending on the values of the parameters. To 
study the stability of the static solution with throat radius clq (if present) we obtain for this metric 
that A is again given by Eq. (|33p . Using that ag > 3M/2, it is straightforward to see that A is 
a positive number, therefore, following Sec. El the static solution is unstable (saddle equilibrium 
point). 

4.3 Schwarzschild-anti de Sitter case 

In the Schwarzschild-anti de Sitter metric the function /(r) is given again by Eq. (|34p . but now 
with a negative A. The event horizon for this metric is placed at 

2/3 



1 - -3^[A|M + ^l + tyAlM 2 



Th = -f — \ i/3 - ( 4 °) 

|A| ( -3^/\A\M + v/l + 9|A|M 2 ' 
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Figure 3: Schwarzschild-anti de Sitter case: number of static solutions for given parameters M, 
A < and A = A + A/(12vr 2 ). Re gion I: no static solutions; region II: one static unstable solution; 
region III: two static solutions, one stable and the other unstable; region IV: one static stable 
solution. 

The horizon radius is a continuous and increasing function of A, with values in the interval 
< < 2M, with r/j(A — > — oo) = and r/ l (A — > 0~) = 2M. The wormhole throat radius ao 
should be greater than r^. The energy density and the pressure at the throat are given by Eqs. 
(137]) and ([38]), with the throat radius a satisfying Eq. ([39]). Using A = A + A/(127r 2 ), which now 
can be positive, zero or negative, we have that Eq. (|39|) takes again the form of Eq. (|28|) with A 
replaced by A. Its solutions are shown in Fig. [2j For the stability analysis, it is not difficult to 
see that A is again given by Eq. (f3"3"|) . Then we have five possible situations, depending on the 
different values of A: 

1. When AM 2 > (547T 2 )" 1 there are no positive real solutions of Eq. (I39p . so we have no static 
solutions. 

2. When AM 2 = (547T 2 ) -1 , Eq. ([39]) has one positive double real root g$ = 3M/2 for which 
A = 0, then following Sec. it is unstable. For the existence of the unstable static solution 
a large enough value of |A| is needed, so that Oq > r^. 

3. When < AM 2 < (547T 2 )- 1 , the solutions of Eq. ([39]) are given by Eqs. (^9ll3ll . with A 
replaced by A (see Fig. E}. Using that a% > 3M/2 and M < ag < 3M/2, it follows from Sec. 
[3] that the static solution with throat radius Oq is unstable (saddle equilibrium point) and 
the solution with radius ag is stable (center). This stable static solution exists if |A| is large 
enough so that a s Q > r^. 

4. When AM 2 = 0, Eq. ([39]) has only one real solution given by = M. The associated 
wormhole solution exists if |A| is large enough so that M > and in this case it is stable 
(center) because A is negative. 
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5. When AM 2 < 0, we have that Eq. (|39p has only one real solution given by 



-\ 2 / 3 



1 - ( / -3vry // 6^i|M + I + 54tc 2 \A\mA 

" -^7i> ( 41 ) 



2itJ6\A\ -3irJ6\A\M + Jl + 54ir 2 \A\M 2 



which is an increasing function of A and lies in the range < aft < M (see Fig. [2]). Then 
we have that A < and the solution is stable. Again, |A| should be large enough to have 
aft > r h . 

The number of static solutions for the different values of the parameters is shown in Fig. 
When AM 2 < (547T 2 ) -1 there are static stable solutions if |A| is large enough so the condition of 
aft > is satisfied, which corresponds to the regions III and IV of Fig. [3j 

4.4 Reissner-Nordstrom case 

The Reissner-Norsdtrom metric represents a charged object with spherical symmetry which has 

/M-l-^ + £ (42) 

where Q is the charge. For \Q\ < M this geometry has an inner and an outer (event) horizon given 
by 

r ± = M ± ^JM 2 - Q 2 , (43) 

if \Q\ = M the two horizons merge into one, and when \Q\ > M there are no horizons and the 
metric represents a naked singularity. When \Q\ < M the throat radius ao should be taken greater 
than r/j = r + so that no horizons are present in M. If \Q\ > M the condition ao > assures that 
the naked singularity is removed. Replacing Eq. (|42f) in Eqs. (]14p and (I15p . we obtain the energy 
density and pressure at the throat: 



q = - v "° ZT , (44) 



y/a 2 - 2Ma + Q 2 

' 1 

(45) 



27rAal 



y/a 2 - 2Ma + Q 2 

Replacing the metric (J32J) in Eq. (fl~3|) . it is not difficult to see that the charge cancels out and the 
throat radius should satisfy the cubic equation (I28p again, with its solutions given by Eqs. (|29H3ip 
and plotted in Fig. [TJ As pointed out in Sec. 14.11 the number of roots of the cubic (|28p are zero, 
one or two depending on if AM 2 is, respectively, greater, equal or smaller than (547r 2 ) -1 . The 
solutions should satisfy that ao > r^, then the number of static solutions will also depend on the 
value of the charge. As is a decreasing function of \Q\, for large values of charge there will be 
two static solutions with radius Oq and aft. Following Sec. El the relevant quantity for the analysis 
of the stability is A, which is again given by Eq. |J55j>. Using that when AM 2 < (547T 2 )- 1 we have 
a s < 3M/2 and aft > 3M/2, it is easy to check that A < for aft and A > for aft, therefore 
aft is stable (center) and aft is unstable (saddle). The special case where AM 2 = (547r 2 ) -1 and 
aft = aft = 3M/2 is also unstable because A = 0. The number of static solutions depends on the 
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Figure 4: Reissner-Nordstrdm case: number of static solutions for given parameters A, M and Q. 
Region I: no static solutions; region II: one static unstable solution; region III: two static solutions, 
one stable and the other unstable. 



values of the parameters A, M and Q. Using that Oq and should be greater than r^, given by 
Eq. ([4*3]) . and defining the functions 

/ , . \ 2/3 

1 + -3irV6AM + iy/1 - Mir 2 AM 2 



and 



« = — rm > ( 46 ) 

2ttVQAM (SttVqAM + Wl - 54tt 2 AM 2 



1 + iV3 - (1 + is/3) (-3ttVQAM + iy/1 - 54ir 2 AM 2 ) 



2/3 



= — ,1/3 '— > ( 47 ) 



it is easy to see that the three possible cases are: 

1. N o static s olutions: when AM 2 > (547T 2 )" 1 , or if a < AM 2 < (54vr 2 ) _1 and \Q\/M < 
\/a(2 — a), where ao ~ 1.583 x 10~ 3 . 



2. One unstable static solution: when < AM 2 < a and \Q\/M < \J/3(2 — 0), or if q < 
AM 2 < (547T 2 )- 1 and yja{2 - a) < \Q\/M < v / /3(2-/3). 

3. Two static solutions, one stable and the other unstable: when < AM 2 < (54-7r 2 )~ 1 and 
\Q\/M > v73(2 -/3). 

The three regions are plotted in Fig. [H For \Q\ > M and ^4M 2 < (547T 2 ) -1 there is always 
one stable static solution (center). Also one stable static solution can be obtained with values 
of \Q\ smaller than M when AM 2 is slightly smaller than (547r 2 )~ 1 . These stable configurations 
correspond to values of the parameters within the region III of Fig. HI A phase portrait of curves 
surrounding a static stable solution is shown in Fig. 
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Figure 5: Reissner-Nordstrom case: phase portrait around a stable static solution. The values of 
the parameters used in the plot are \Q\/M = 0.95 and AM 2 = 1.85 x 10 -3 and the radius of the 
throat is ag/M = 1.406. The perturbation e(r) is defined by e(r) = a(r)/ag — 1 and u(t) = e(r) is 
its time derivative. The curves have different initial conditions and rotate clockwise. 

5 Conclusions 

We have constructed spherical thin-shell wormholes supported by exotic matter fulfilling the 
Chaplygin gas equation of state. Such kind of exotic matter has been recently considered of particu- 
lar interest in cosmology as it provides a possible explanation for the observed accelerated expansion 
of the Universe. For the wormhole construction we have applied the usual cut and paste procedure 
at a radius greater than the event horizon (if it exists) of each metric. Then, by considering the 
throat radius as a function of time we have obtained a general equation of motion for the Chaplygin 
gas shell. We have addressed the issue of stability of static configurations under perturbations pre- 
serving the symmetry. The procedure developed has been applied to wormholes constructed from 
Schwarzschild, Schwarzschild-de Sitter, Schwarzschild-anti de Sitter and Reissner-Nordstrom ge- 
ometries. In the pure Schwarzschild case we have found that no stable static configurations exist. A 
similar result has been obtained for the case of Schwarzschild-de Sitter metric (positive cosmological 
constant). For the Schwarzschild-anti de Sitter geometry, we have found that the existence of static 
solutions requires that the mass M, the negative cosmological constant A and the positive constant 
A characterizing the Chaplygin fluid should satisfy AM 2 < (547T 2 )" 1 , with A = A + A/(127r 2 ), and 
|A| great enough to yield a small horizon radius in the original manifold. If these conditions are 
verified, for each combination of the parameters, when < AM 2 < (54-7r 2 )~ 1 there is one stable 
solution with throat radius ciq in the range M <a < 3M/2, and for AM 2 < there is one stable 
configuration with ao < M. In the Reissner-Nordstrom case the existence of static solutions with 
charge Q requires AM 2 < (547T 2 ) -1 ; we have found that if this condition is satisfied, then one 
stable configuration always exists if \Q\/M > 1. When \Q\/M is slightly smaller than one, there is 
also a stable static solution if AM 2 is close to (54-7T 2 ) -1 . In this work, then, we have shown that 
if over-densities in the Chaplygin cosmological fluid had taken place, stable static configurations 
which represent traversable wormholes would be possible. 
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